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1. Introduction 



In this paper we discuss a solution for the high energy scattering amplitude in the case of high density QCD. We 
would like to go beyond the mean field approximation, where the solution has been discussed and well understood 
both analytically and numerically (see [0, 0, & 0, plj^, Sl) , and to approach the summation of the so called Pomeron 
loops that should be taken into account |Q, |o| [ll|, p^, LSI. The problem of taking into account the Pomeron loops 
can be reduced to the BFKL Pomeron calculus |1J, ^ |l6| and/or reduced to the solution of the statistical physics 
problem, i.e. the Langevin equation and directed percolation [0 ^ 

At the moment the problem of the high energy amplitude in QCD is well understood and can be solved only in 
the kinematic region given by 



1 < ln(l/a|) < asF < — 

as 



(1.1) 
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To go beyond this region we need to know the corrections of the order a| to the BFKL kernel as well as the corrections 



to the vertices of the Pomeron interaction. Eq. (1.1) will play a significant role in our discussion and we wish to 
present a more detailed derivation. It is well known that the exchange of one BFKL Pomeron leads to a contribution 
which can be written as 

A (one Pomeron exchange) cc a'^e^^ (1-2) 

where Y = ln(s/so) ( s = W'^ and W is the energy of the scattering in the cm. frame ). The intercept of the BFKL 
Pomeron can be represented in the form 

A = as XLO BFKL (7 = 1/2) + a| XNLO BFKL (7 = 1/2) (1-3) 

where x is the Mellin transform of the BFKL kernel (see more details in ). 
The contribution from the exchange of two BFKL Pomerons is proportional to 

A (two Pomeron exchange) oc (a^e'^^)^ (1-4) 

Comparing Eq. (|l|) to Eq. ^) one can see three different kinematic regions: 

1. 2 ln(l/as) > Ay > 1 ; 

In this kinematic region we need to take into account the BFKL Pomeron at leading order since a^Y 1 
and neglect the contribution from multi Pomeron exchanges. 

2. 1/as > Ay > 2 ln(l/a5) ; 

We need to take into account multi Pomeron exchange but we can still restrict ourselves to the BFKL kernel 
in the leading order approximation. 

3. Ar > 1/as; 

In this region we have to calculate the next-to-leading BFKL corrections to the BFKL kernel as well as the 
corrections to the Pomeron vertices. In addition, it is not clear whether or not we can rely on the BFKL 
Pomeron calculus in this region |21| . 



Concentrating on the kinematic region of Eq. (1.1) we develop a strategy which consists of three steps. First, 
we show in this paper that the BFKL Pomeron interaction in this kinematic region can be reduced to the exchange 
of non-interacting Pomerons if we renormalize the low energy amplitude of the interaction of 'wee' partons with the 
target. 

Using this observation we propose an improved Mueller-Patel-Salam-Iancu method (MPSI) for summation of 
Pomeron loop diagrams. This method is actually the i-channcl unitarity constraints re-adjusted in a convenient form 
for use in the dipole approach of QCD p3| . 

Finally we propose that the answer which we obtain, is the real solution to our problem at ultra high energy. 
We will justify and show that this is a reasonable hypothesis. 

This paper is organised in the following way. In the next section we discuss the simple case of the BFKL Pomeron 
calculus in zero transverse dimensions. Being the simplest model for the Pomeron interaction this approach allows 
us to discuss our main ideas and suggestions without complex calculations. Actually, the main content of this section 
has been discussed in p^ , but for completeness we present this model in our paper. We hope that the reader will 
find this instructive in later sections. 
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Figure 1: The renormalization procedure in tlie case of tlie simplest 'fan' diagram. 



In section 3 we argue that the main properties of the BFKL Pomeron calculus in zero transverse dimension are 
inherent for the BFKL Pomeron calculus in QCD. In particular, we can consider the scattering amplitude in the 
kinematic region of Eq. (1.1) as the exchange of BFKL Pomerons neglecting their mutual interactions. 

Section 4 is devoted to the calculation of the scattering amplitude in the model for the BFKL kernel which has 
been developed in Ref. [ p5| . In this section we demonstrate how the MPSI approach works and we obtain a formula 
for the scattering amplitude. We advocate that the resulting scattering amplitude satisfies the unitarity constraints 
both in the s and t channels and could be a good candidate for the answer outside the kinematic region given by 
Eq. (0). 

In our conclusion wc summaries our results and compare them with the approaches that we currently have on 
the market. 



2. The BFKL calculus in zero transverse dimension: non-interacting Pomerons and 
improved Mueller-Patel-Salam-Iancu approach 

In this section we analyse the Pomeron diagrams in the BFKL Pomeron calculus in zero transverse dimensions (toy 
model of ^). Recently, this model has become rather popular (see and |^ where the most important 

aspects of this model have been discussed and solved) 

Firstly we consider the simplest 'fan' diagram of Fig. and Fig. |. It can be calculated in an obvious way, using 
the explicit expression for the Pomeron Green function, namely, G{Y — Y') — exp {A {¥ — Y')). Indeed, for the 
diagrams of Fig. |^ we have 

A{Fig.^ = 7G(r-0) - Aj^ / dY' G{Y - Y') G^{Y' - 0) (2.1) 

= je^^ - A72 rdy'e^(^+^') (2.2) 
Jo 

= le^^ - A7^ (le^^^ - ^e^^) 

= _ ^2 g2Ay ^ _^ ^2) gAl- ^ _ ^2 g2Ay ^ ^A F 

where A is the Pomeron intercept which is equal to the triple Pomeron vertex r(l — )■ 2) in this oversimplified model, 
and 7 is the amplitude of the Pomeron interaction with the target. 

For the diagrams of the second order, given by Fig. H, we have to integrate over the two rapidity variables yi 
and 2/2- The result is 
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Figure 2: The renormalization procedure in the case of the 'fan' diagram of the second order. 



A{Ftg.§ = 

= 2A^j^ r dyi r dy2G{Y-y,)G{yi-0)G{y2~0)G^{y2-0) (2.3) 
Jo Jo 

= 2A'j^ J\y,£' dy,e^(^+y^+y^^ = 2A^3(^^,3Ar _ ^^2Ay ^ _^^Ay 

As one can see, the integration over Y' in Fig. |l| and over yi and y2 in Fig. reduces these two diagrams to 
the sum over diagrams which describes two contributions: the exchange of two non-interacting Pomerons and the 
exchange of one Pomeron with the renormahzed vertices: ^^j^"^ —7 + 7^ ^nd 7^^'' =7 + 7^+7"^ 

Adding the two contributions of these diagrams gives 

A{Ftg.^ + A{Ftg.§ = (2.4) 
= ~7'e3Ay „ (^2^273)6^^^ + + ^2 ^ 73)e3'^^ = 7' e^^^ - {j^i^f e^^"" + 7^6^^ 

Therefore, one can see that the scattering amphtude can be rewritten as the exchange of Pomerons without interaction 
between them but with the renormahzed Pomeron - particle vertex. In the dipole model this vertex is the amplitude 
of the two dipole interaction in the Born approximation of perturbative QCD. 

These two examples illustrate our main idea: the BFKL Pomeron calculus in zero transverse dimensions can be 
viewed as the theory of free, non-interacting Pomerons whose interaction with the target has to be renormahzed. 
The master equation for the scattering amplitude in the mean field approximation can be easily rewritten (see Ref. 
[p3) in the form 

oY djR 

with 

IB. = (2.6) 

1 - 7 



Thus we have shown in Eq. ( |2.lD and Eq. ( |2.3| ) how Eq. (2.6) has started to build up in the perturbative expansion. 

The general solution of Eq. ( |2.5| ) is the system of non - interacting Pomerons and the scattering amplitude can 
be found in the form 

00 

No {jr\Y) = ^ (-1)" C„ 7fl G"(y - 0) (2.7) 

n=l 

where the coefficients C„ could be found from the initial conditions, namely, from the expression for the low energy 
amplitude. In particular, the initial condition 

NohnlY^O) = 7 = iR/a + lR) (2.8) 
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generates C„ = 1 and the solution is 



NoilRlY) = ^ (2.9) 

The initial condition of Eq. ( |2.8| ) has very simple physics behind it that has been discussed in Ref . ||2^ . 

For the analysis of the enhanced diagrams we start from the first diagram of Fig. ^, which can be written as 
follows 

A{Fig.§ = (2.10) 

-A272 f dy, r dy2G{Y^y,)G^y,~y2)G{y2-0) 



Jo 



= 1 dy, J' dy2e^^^+y^-y^'> 

where r(2^1) = A72 (see Rcf.|li). 



Adding Eq. (2.10) to the exchange of one Pomeron we obtain the following 

One Pomeron exchange + A (i^ig. I) = jRe^"^ - j'^e'^^^ (2.11) 

with 

IR = 7^^^ = 7 + 7^ Afi = A + 7A; (2.12) 

Therefore, the Pomeron loops can be either large ( where their size in rapidity is of the order of Y) and they can be 
considered as un-enhanced diagrams; or small (where their sizes are of the order 1/A) and they can be treated as 
the renormalisation of the Pomeron intercept. 

In QCD A oc as while 7 oc a^. Therefore, the renormalisation of the Pomeron intercept A is proportional 
to a|. We can neglect this contribution since (i) there a lot of a| corrections to the kernel of the BFKL equation 
that are much larger than this contribution; and (ii) in the region of y ^ l/o^s ' where we can trust our Pomeron 
calculus (see introduction) (A/j ~ A)Y 1. 

Concluding this analysis we can claim that the BFKL Pomeron calculus in zero dimension is a theory of non- 
interacting Pomerons with renormalised vertices of the Pomeron-particle interaction. In the dipole language, it means 
that we have a system non-interacting Pomerons with a specific hypothesis on the amplitude of the dipole interactions 
at low energy. For the problem that we are solving here, namely, when we have one bare Pomeron at low energy, 
this amplitude is determined by Eq. (2 



For such a system we can calculate the scattering amplitude using a method suggested by Mueller, Patel, Salam 



and lancu and developed in a number of papers (see Refs.|23, 3™ B2l pu pa, pq| and references therein). In this method 



the scattering amplitude is calculated using the t-channel unitarity constraints which is written in the following way 
(assuming that the amplitudes at high energy are purely imaginary, N — Im A) : 

Ni[...]\Y) = N{[...]\Y-Y';P^nP)(^N{[...]\Y';P^nP) (2.13) 

where stands for all necessary integrations while [. . . ] describes all quantum numbers (dipole sizes and so on ). 
The correct implementation of Eq. ( ^.13 ) leads for our case to the following formula (see also Rcfs. |^, ^ pof ) 



^^0 



MPSI 



m - 1 -»p|-7''^^^|»''^^(h!,''l>--r')iV''-((7g'|V-')l,™.,™.„ (2.14) 
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where N'^'^^^ {Y^^b) is given by Eq. (2.£)(see also Eq. (2.5)) in the mean field approximation and 7^" oc is 



(X a| 

the scattering amplitude at low energies which is described by the Born approximation in perturbative QCD. The 



difi^erence between Eq. (2.14) and the original MPSI approach is the fact that this equation does not depend on the 
value of Y' and, because of this, we do not need to choose Y' = Y/2 for the best accuracy. 



Substituting Eq. (2.9) in Eq. (2.14) we obtain 



MPSI (^,BA 



= 1 - exp 



1 



BA„T(1^2)Y 



1 



7 



BApT{1^2)Y 



r 0, 



1 



(2.15) 



r (0,a;) is the incomplete gamma function (see formulae 8.350 - 8.359 in Ref. Gol 



We claim that Eq. ( 2.15 ) is the solution to our problem. One can easily see that Nq (7IF) ^ 1 at high energies 
in contrast to the exact solution (see |^). The exact solution leads to the amplitude that vanishes at high energy 
(see Refs.pgj, ^). As has been mentioned the solution depends crucially on the initial condition for the scattering 
amplitude at low energies. For Eq. ( 2.15| ) this amplitude is equal to 



^MPSI 



(7|F = 0) = £(-l)"+in!(7 



BA\ 



(2.16) 



n=l 



with 7 



BA 



Eq. (2.14) can be rewritten in a more convenient form using the Cauchy formula for the derivatives, namely. 



n yMFA 



ilR\Y) 



97S 



2 7ri 



(2.17) 



The contour C in Eq. (2.17) is a circle with a small radius around 7^ — 0. However, since the function Z does not 
grow at large 7^; for n < 1 we can close our contour C on the singularities of the function Z. We will call this new 
contour C^. 



N, 



MPSI 



{Y) = l - exp<^ -7 



BA 



d d 



= 1 



^ Hn' Hp' 

{-l^^T . . 1 



^MFA 1^(^(1) |y _ Y') N^'^^ ((7^'^') 





zMFAuM) 



h'n'lY-Y') 



cl 



J (2) 

dip 



7MFAf„i^) 



1 — exp 



ryBA^y 



.(1) ,(1) / -(l)-(2)- 

7i? 7_R p [q Tr Ip 



ryBA^y 



jBA^y 



X Z 



MFA ~. 



(7«) {^f) 



Here we introduce new variables 7]^ 

r,MFA 



7^'^ exp {y - y) and 7^,'^ 



71,^^ - 



Closing the integration on the poles 7]^ 



-1 and 7^^^ — 



ZMFA = 



7]j,^ exp {y'). In these new variables 
1 



1 + 7i^^ 



(2.18) 



(2.19) 



-1 we obtain the formula of Eq. (2.1 



3. Calculation of the simplest diagrams in the BFKL Pomeron calculus 



The main goal of this paper is to show that in the general case of the BFKL Pomeron calculus we have the same 
situation as in the simple model. We start from the analysis of the simple diagrams of the BFKL Pomeron interaction. 
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Figure 3: The renormalization procedure in tlie case of tlie simplest eniianced diagram. 



3.1 The functional integral formulation of the BFKL Pomeron calculus 



The most economic form of the BFKL Pomeron calculus was suggested and developed by Braun in Ref.|16 where 
he formulated this theory as the functional integral 



= j D^D^+e^ with^ = Sq + Si + Si 



(3.1) 



where describes free Pomerons, Si corresponds to their mutual interaction while Se relates to the interaction 
with external sources (target and projectile). 

They have the form 



(3.2) 



Si 



2nai 



dY 



d x^-y^ d x_'2 d X 



Li2 ±.23 ±31 



{(Li,2<f(xi,a;2;r)) <i>+(^2,3l3;n$+fe,2Ii;r) + h.c.} ; (3.3) 



Se = - dY d'^Xid'^x2{^{xi,X2;Y)Tpr{xi,X2;Y) + (xi,X2;Y) Ttar{xi,x2;Y)} 



(3.4) 



It should be stressed that Eq. (3.4) describes the local interaction both in rapidity and in coordinates and Tpr (rtar) 
stands for the projectile and target, respectively. 

The Pomeron Green function has the form Q 

G {Xi, X2',Y\Xi, X2',Y ) = 

9{Y-Y') J2 / di^ e"("''^)(^~^)(zy2^— )A(n,:/)i;*(a;'i,£'2;a;o|i/)£;(xi,a;£ok) (3.5) 

n— — oo 

where the vertices E are given by 



Eix-i^,X2;xo\n;iy) = {-If 



Xl2 
XlO X20 



^12 



(3.6) 



; /i = (1 + n)/2 + iv and h = I - h* = [l - n) /2 + ii 



^Xi^x and Xi^y are x- and y-components of the two dimensional vector x^ 
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The energy levels are the BFKL eigenvalues 



\n\ + l 



(3.7) 



uj{n, v) — as \y4>{^) — Reip 

with ipiz) = dhiT{z)/dz and T(z) is the Euler gamma fmiction. 
The operator L12 in Eq. ( ^T| ) is defined as 

Li,2 = 3Ii2P?P2 with £2 = -v^ (3.8) 

with E{xi,X2;XjoW) ^^'^ x{n,u} being its eigenfunctions and eigenvalues respectively. It is easy to check that 

1 



with 



Li,2E{xi,X2;Xf)\n;iy) 
X{n, v) 



X{n, v) 
1 



E{x^,X2]XQ\n;u) 



(«+i)^ 



91 r (n— 1)2 



For our further discussions we rewrite Eq. (3.6) in a more convenient form as 



2^12 



''12 



,(^12-^)(Z12 + ^) 

where the new variable Z12 is defined by 

Xl + X2 

Z12 = — :z xo 



, 1^12 2 ^-^12 ~^ 2 ) 



(3.9) 
(3.10) 



(3.11) 



(3.12) 



3.2 Triple Pomeron vertex in momentum representation 

In this section we calculate the triple BFKL Pomeron vertex which is the key ingredient of all enhanced diagrams of 
the type illustrated in Fig. |[ We find it easier to do all the calculations in the momentum representation. We define 
the Fourier transform of the vertex function in Eq. (3.6) as 



g(fc, g, n, v) 



(Px 



<fze'^-'^^^ E{x^,X2]x^\n]v) 



(3.13) 



with and z = {x-^+X2)/2 — XQ = b — Xj^. As it follows from Eq. ( t3.13| ), the momentum k is defined 

as the momentum conjugate to the dipole size. It will be shown later that the physical meaning of the momentum q 
is a momentum transferred along the Pomeron. The explicit expression for the function g{k, q,n, v) is found in the 
Appendix A. 



Using the definition of Eq. (3.13) we want to express the free Pomeron Green function defined in Eq. ( p.5| ) in 
terms of the function g(k, q, n, i^). To do this we find an inverse Fourier transform of g{k, q, n, v) and substitute it 
into Eq. (3.5). The vertex function E in terms of g{k, q, n, v) reads 

(Pk dPq 



E{x,^X2;x^\n-v) = x,2 J ^ 



e '--12 e" 



''--g{k,q,n,i^) 



(3.14) 
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With the help of Eq. we rewrite the free Pomeron Green function in Eq. ( |3.5| ) as 

(X-i, X2',Y\Xi, X2',Y ) = 

oo „ 2 



(3.15) 



1_ p+tk'-x[^ g+V-£'i2 „*(k'.n' 



(27r)2 (27r)2 



lik, q, n, v) 



n— — oo 



„'2 



(27r)2 (2^)2 ^ 9 {K,q,n,y) x^2 



(27r)2 (27r)2 



g --12 e e ■ g(fc, g, n, i/) 



The integration over brings in a delta function (27r)2(5(2) _ q'y After the integration over the momentum q' we 
end up with the expression for free Pomeron Green function in terms of g{k, q, n, v) as follows. 



G {Xi, X2',Y\j/i, i/2tY') — 

oo „ 2 

9{Y -Y') J2 / e"^"''^)^^-^') (i.2 + !L) A(n, i.) (2^)2^2 _ ^) x 



(3.16) 



n—~oo 



±12 



(27r)2 (2^)2 ^ - 9 {!£,q,n,,^)x,2 



Pk (fa 

(2^)2 (27r)2 9(K,q,n,iy)- 



oo „ 2 

eiY-Y') ^ / d7.e"("''^)(^-^')(z.2 + !L)A(7^,7.) X 



±12 ±12 



n— — OQ 
Pk' dPk Pq 



9*ik,q,n,i^)g{k,q,n,iy) 



(27r)2 (27r)2 (27r)2 

The expression in Eq. ( |3.16| ) clarifies the meaning of the momentum q defined in Eq. ( |3.13| ) as the momentum 
conjugate to the difference in the impact parameters of the interacting dipoles or, alternatively, as the momentum 
transferred along the Pomeron. 

Using the same standard procedure we want to calculate the triple Pomeron vertex defined in the BFKL Pomeron 
Calculus by Eq. (3.3). It is better to write the vertex in terms of free Pomeron Green functions as 

dY 



2 7ra| 



d^x.,d^^2d^x, ^^^^^ (ai;,a^; yi|ai,£2; Y))G ix2,x,; Fla^.a^i; ^2) G ix,,x,;Y\x!„x:,;Ys) (3.17) 



^12^23^31 



and then use the representation of G (x'j^, 2:2; K) in terms of the vertex function E defined in Eq. (3.5). 

Thus one can see that the calculation of the triple Pomeron vertex in terms of the functions g{k, q, n, v) is reduced 
to the calculation of 



2 nag f d'^x^d^X2d x 



Nr. 



£^12^23^31 



■(-^12^^(^l,£2;^0ll"'i;'^l))^*(£2'23;^02l"2;i'2)£^*(^3,^i;£03l'^3;i^3) 



(3.18) 



then being properly convoluted with other necessary vertex functions E as defined in Eq. (3.5) to reproduce Eq. (3.17). 



We use the fact that the vertex function E are eigenfunctions of the operator L12 (see Eq. (|3.9D) to recast Eq. (3.18) 



as 



£12^23^31 A(71l,^l) 



£^(2l,22;20ll"-i;^l)^*fe2'23;£02l"2;i^2)^^*(23,£i;203l'^3;J^3) (3.19) 
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Following the standard procedure we used above in finding the Pomeron Green function in terms of functions 
g(k, q, n, v) we insert the inverse Fourier transform of Eq. ( |3.14 ) into Eq. ( 3.19| ) and then integrate over the coordinates 
, Xo and follows. 



2 Tral 

Nr. 



!i23 



L31 



{2nf 


(27r)2 








(27r)2 






(27r)2 


(27r)2 



12 -£12 e ^-ii^ (7(fci2,9i,«i,i^i) 



(3.20) 



2 7ra| 
Nr 



<fx^(Px2(Px^ f (fk^2 '^'^1 



(27r)2 


(27r)2 






(27r)2 


(27r)2 




d\ 


(2^)2 


(27r)2 



■11+^2 



+m -a; 
e -1 - 



■01 



X 



^ g+»fc,3.fe-x3) g+»i,-=-^ e-*^2-2o2 g*{k23,q^,n2,iy2) X 



. ^2+^3 



e-'^3-^ g-»23-io3 g*{k^^,q^,n3,P3) 



The integration over coordinates in Eq. (3.21) gives 

/ d'x, =^ (2^)2^(2) ik,2 + \q^ + - \q^) 



Jd\2 
Jd^Xs 



{2n)H^^^){-k,2 + \q,-k2,-\q^) 
(2^)'^('Hfc23-k2-^31-k3) 



Now we perform the integration over fcjs: .^31 ^-nd q which results into 



(3.21) 



(3.22) 



/ d^k23 



^31 = k3 " -^12 - 111 '^'^^ (2vr)25(2) (ig^ - - kl + k2 " ^2 

^23 = k3 " -^12 - kl + k2 ^'^^ (27r)2,5(2)(g^ - - 



93 ^^1-^2 



(3.23) 



where we used the identity / dx f{x) d{x — a) S{x — 6) = /(a) S{a — b). The result of Eq. (3.23) can be summarized 
as 

1 



i31 



-k 



12 2^2 



-^23 ~ ^12 + 2-2 



(3.24) 



^3 =ll-l2 



Finally, we are in position to write the triple Pomeron vertex in terms of functions g{k, g, n, v). To do this we remove 
the performed integration over x^, X2, x^, k^i, A:23 and q^ in Eq. (3.21) as was shown above by replacing k^^, ^23 and 
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q with the expressions given in Eq. (|^), as well as, multiply Eq. ( |3.2lD by {2Trf from (5-functions. Thus Eq. ( ^.2lD 
after the outlined integrations reads 

d2fc d^q d?q^ 



(2^ (2^^ (2^ e-i.-.-'i. .o.-.,.-..)-.03 X r (3.25) 
where 

fi \ lirai 1 ,^,,11 \ */ , 1 N 



The expression in Eq. (3.25) represents the triple Pomeron vertex in terms of functions g{k^ q, n, v). As it was shown 
on the example of the free Pomeron propagator, the triple Pomeron vertex should be properly convoluted with all 
the necessary vertex functions E to give the amplitude of a dipole being scattered off two dipoles. The only missing 
thing is the explicit expression for g{k^ q, n, v) which is calculated in the Appendix A and should be plugged into 
Eq. ( 3.16| ) and Eq. ( 3.25| ) for obtaining the final expression. In the region of 4fc2 ^ q^ the ftmction g{k,q-,^,v) is 



independent of q^ and has the form (see Appendix A) 

9{kqM^C{v){er"^+^^ (3.26) 

where the constant C(i^) is given by 

—IV 1 ^(1/2 + ii/)I ^^^) 

3.3 The simplest 'fan' diagram 

In this section we calculate the set of diagrams shown in Fig. |l|, namely, single Pomeron exchange and the first 'fan' 
diagram for the triple Pomeron interaction. 

3.3.1 The single Pomeron exchange 

According to our definition of the function g{k,q,n,iy) given in Eq. ( 3.13| ) the momenta k and q denote momenta 



conjugate to the size and the impact parameter of the dipole respectively (see also the remark after Eq. ( 3.16 )). Using 
this notation we readily write the expression for the dipole with rapidity Y and the transverse momentum k being 
scattered off the dipole with rapidity Y' and the transverse momentum fcg due to the exchange of one Pomeron. 

A{1P) ^ I g*{k,q,0,i') v^ X{0,v) g{ko,q,0,iy)e'^^'''^^^-^"> 

,u.o{'^)(Y-Y') ^ 327r6p(fco; k\Y - Y') (3.28) 







[ '^'^ k ko 







where u!o{v) ~ 4 In 2 — 14C(2)i/2 is the expansion of the BFKL eigenvalue in the vicinity oi v = bringing the main 
contribution at high energies, and P{ko; k\Y — Y') is well known expression for single Pomeron exchange given by 

Piko;k\Y-Y') = l f^-Lf^Y^-^C^^^^-^'^ (3.29) 
TT J 2n k ko \ k J 



The factor of 32 merely reflects our normalization of functions g{k, q, 0, v) defined in Eq. ( |3.13| ) . From Eq. ( |3.2£ ) one 



can see that the anomalous dimension of the gluon structure function {xG{x, Q"^) oc {Q'^)'') is equal to 7 — ^ ~ iv 
since P{kQ;k\Y -Y') cx {l/k'^)xG{x = eyi-piY -Y'),e /kl). 
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Figure 4: The first 'fan' diagram. Wave lines denote tlie BFKL Pomerons. 



3.3.2 The first 'fan' diagram 

The first fan diagram, for interaction of one dipole with two dipoles, is shown in Fig. ^ Using the representation 
for triple Pomeron vertex in momentum space calculated in Eq. ( 3.25| ) we can write the amplitude corresponding to 
Fig. for 71 = as 



AiP^2P;FigM) = 



2T:a 



2 fY 



N, 



c Jo 



dY' 



, f dPUi ^'22 



(27r)2 (27r)2 (27r)2 

.2 ,.2 ,.3 



dv^ dv2 di^3 {2iif6^^\q^ - (2^)2,5(2) _ X 



-1 -2 A(0, j^i) 

g.(.3)(y'-o) 0,7^1 ) 5(fc'i,9;,0,i.i) g*(-fc'i + - ^g^,?;, 0,1^2) 5(^2,92, 0,^/2) x 

5*(-fc'i - ^92' 2i -92' ^^,^^3) 9{k3,q^, 0,1^3) 



where ^2 = ^h!i + 59']^ ~ 592' —3 ^ ^—'1 ~ 5''2 '^3 ^ '^1 ^ '^2 shown in the derivation of Eq. (3.25). 

For the case of small momenta transferred along the Pomeron we approximate the vertex function by expression 
given in Eq. (3.2(:). It is worth to mention that the constants C(i^) defined in Eq. (3.27) have a property 



for 



(3.31) 



Thus Eq. (|^ reads 
-AiP^2P;Ftg.g) 



2T:a 



2 pY 



dY' 



, f d^k!, d^^ 



(27r)2 (27r)2 (27r)^ 



dl^l dv2 dv3 (27r)2^(2)(^^ _ ^, ) (2^)2^(2) (^_^ _ x (3.32) 



(2^)25(2)(2^^_g;+2^) (7r4)3A(0,i/2) A(0,i/3) e-(''^)(^-^') e-(-^)(^'-o) e-t''-^)^^'-") x 



|fciri-2.., |fc;ri+2-^ |-fc; + ,;/2-9;^/2r^-'^-|fc2r'+'^'^^ 



In Eq. (3.32) momenta q^ and q^ are arbitrary. However, in the mean field approximation (MFA) the momenta 
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transferred along the Pomeron are much smaller than momenta conjugate to dipole sizes ^ |fc| ^ which allows us 
to simply Eq. ( 3.32| ) using | — fc'i + — ~ l^i | and | — fc'^ — g^/^l ~ l^'i I • It is easy to see that the integration 



over 'k!_^ leads to a delta function &(i + 2v\ — Iv^ — 21/3), which is the conservation of the anomalous dimensions 
(7 = 1/2 — iv^ in the triple Pomeron vertex (5(1 + 71 — 72 — 73) (see Ref. [ p^ where this (5-function was found for the 
forward scattering amplitude). However we need to be careful since in Eq. ( [3.32 ) we implicitly used the condition 



that k\ < k'l < kf. For fc\ > kf the second term in Eq. (A-21) gives the main contribution. Therefore we integrate 
over all the kinematic region from k'l = to k'l = 00 but subtract the regions k'l > kf and replace them by correct 
integral. 

For example for k'l > k\ we obtain the following integral 

where I = ln{kf/k'^). The result of the integration leads to 

g(-l-iy2-ii'3)' [ }. _ }. ) (3 34) 

y-l + z2j/i - i2i/2 - 2^3 -1 - i2i/i - i2i/2 - 2^3/ ^' ' 

It is easy to see that the integral over is equal to zero. The behaviour of to^i'i) at large values of vi we regularise 
by adding a small ie term for {uj^iyi) —> lo^vi + ie). Lo{h'i + ie) '^^^"^^°°> — oo providing a good convergence of 
the integral on the large circle. 

Therefore, finally we reproduce the (5(7i — 72 — 73) contribution. 

After further integration over g^, and dvi in Eq. ( |3.32| ) we obtain 



A{P^2P;FtgM) = 



dY' / di.2 ^ {2^f^^^Hg,-g,+g^) (vr^)'A(0,t.2) A(0,z.3) x 



2TTdi 



-1 + 2(1/2 



1^3 I 



-l+2ii.3 



2TTdi 

N, 



dv^dv:, - (2^)'<5(2)(g^-g^ + g^) 



[n^fmu^) A(0,z.3) 



X (3.35) 



l-2ij/2-l-2iy3 1^2!""^'''^*''^ \k 



( 



I + 21W3 / u{^i/2 + V2 + V3){Y-0) 



)){Y-0)^ 



As was already shown in the toy model in Eq. (2T) the integration over rapidity brings in two terms in the last line 
of Eq. ( 3.35| ). The first term corresponds to single Pomeron exchange and the second term corresponds to double 
Pomeron exchange of two non-interacting Pomerons. To see this we rewrite the two terms in a more convenient form. 
The first term in the last line of Eq. ( 3.35|) reads 



27ra2 1 



Nr_ 2 



-(27r)2<S(2)(q^-2^+^^ 



di'2 dn 



{TTy\{0,V2) A(0,z/2 + M~'^2) 
uj{^i) ~ uj{v2) - uj{i/2 + ^1- V2) 



(3.36) 



|^j-l-2iAi l/fcjP"^^^*"^ 1^^1-1+2^-1-211-2 gC^(^)(y-0) 



^To see this It is enough to recall that we can safely use MFA for the nucleus target where fei k2 
where R and Ra are nucleon and nucleus radii and Ra 2> R- 



l/R but gi 92 ~ 
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where fi = —i/2 + 1'2 + 1'3. At high energies the leading contribution conies from smaU values of fi and we can simplify 
Eq. ( p6| ) as 



2nal 1 
^2 



(7r4)3A(0,t/2) A(0,-i^2 + V2) 1 fk 



1 



fcifcs V fci 



2TTal 1 

^vr2 



(3.37) 



where P{k3;ki\Y — 0) is a single Pomeron exchange defined in Eq. ( 3.29 ) and f{k2,k^) is a function of fc2 and k^ 
which has to be calculated and is given by 



/(fe, kz 



(7r4)3A(0,i/2) A(0,-i/2 + i/2) 1 A3 



2ii/2 



(3.38) 



The second term of Eq. (3.35) represents two non-interacting Pomeron exchange. To see this we recast it in the 
form of 



27ra| 
1 



fclfc2 V^l 



1 



dy2dv^-{2TT)H^^\q^-q^+q^) 



(7r4)3A(0,z.2) A(0,i.3) 



u{~i/2 + V2 + V'i) - UJ{V2) - IjjiV'i) 



2iV2 



1 



kik^i \ki 



2lV3 



(3.39) 
(3.40) 



At high energies the main contribution comes from small values of 1^2 and ^3 and thus we can simplify Eq. ( p.39| ) as 
follows 



'"^^^2.)V^)(,3-,^+,^) (4^^ 



Nc 2 



Lu{-i/2)- Lj{0)- Uj{0) 



(27r^)^P(fc2, ki\Y - 0) Pik3,ki\Y - 0) (3.41) 



Comparing Eq. ( |3.35D with Eq. (|2T|) using the last result we conclude that in QCD the 'fan' diagrams can be 
viewed as the contribution from the exchange of two non-interacting BFKL Pomerons and a renormalisation of the 
single Pomeron contribution. 

However, as it was noticed by Hatta and Mueller ji^, there is another saddle point in Eq. ( 3.35 ) where the 
denominator lj{v2 + 1^3 — i/2) — uj{v2) — ^{^3) is close to zero. Indeed, the equation for this saddle point ( jsp = 
1/2 — ivsp^i) looks as follows 



2uo'{^Sp)Y - 



2{U.'{2{^SP-1)-^'{1SP) , 1 .^2/fc2^ _ 

c.(2 7sP-l)-2c.(7sp) +^"(^i/^-3^ - ^ 



(3.42) 



In Eq. (3.42) we assume that the main contribution is dominated by vi — V2 One can see that the position of this 
saddle point is very close to the solution of the equation 

a;(2 7o-l)-2c^(7o) = (3.43) 

Denoting 7sp — 70 by ^7 one can see from Eq. (3.42) that for the first term in Eq. (|3.35 ) we have 

1 



^7 — 



2^(7o)y + ln(fc2/fc|) ^" 



and for the second one 



(57 



2w(27o-l)r ln(fcf/fc|) 



< 70 



(3.44) 
(3.45) 
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Figure 5: Solution to Eq. (3.43) for ti;'s, given by the Figure 6: Solution to Eq. ( |3.43| ) for uojyert given by the 
BFKL equations. BFKL equation and uJaat calculated from Eq. (3.47). 



The sum of these two terms leads to the contribution with energy dependence 



{P^^P-Fig. 



,(w(7o+'i<5i^)+w(7a-i(5i/)) {Y -0) 



It turns out (see Fig. ||) that the resulting intercept 20; (70) is larger than the intercept for the two Pomeron 
exchange (2cl'(7 = 1/2)) and, therefore, this contribution is the largest among the three. The appearance of such new 
singularities in the angular momentum plane is a deadly blow to the entire approach based on the BFKL Pomeron. 
Indeed, since the new singularity (double pole) is located to the right of the singularity generated by the exchange 
of two BFKL Pomerons, we need firstly to sum over all such singularities, to obtain the resulting Green function 
and only after doing this can we build the Reggeon calculus based on this Green function. However, the situation 
changes crucially if we take into account the fact that in Eq. ( |3.43 ) ^(70) and 1^(270 — 1) are actually in quite different 
kinematic regions. a;(27 — 1) enters with a small value of the argument, namely 270 — 1 < 7cr and, therefore, can 
be calculated using the BFKL equation of Eq. (3/7). We recall that 7cr can be found from the following equation 
is 011 

X {lev) dx (7cr) 



(3.46) 



However, we have a different situation for a;(7o). Indeed, 70 turns out to be larger than ^cr (see Fig. ^ and recall 
that 7cr = 0.37). Therefore, 0^(70) describes the behaviour of the scattering amplitude in the saturation region, and, 
therefore, cannot be calculated using Eq. (3/7) (see also Ref.|0], where the same conclusion has been derived from 
slightly different considerations) . As it was found in Ref. for uj in the saturation region we have 



1 - 7cr 



(1-7) 



(3.47) 



at least for 7 > jcr but close to 7cr- 

As one can see from Fig. I Eq. ( |3l3| ) which can be written as 



2Wsai(7o) = t^pert(27o - 1) 



(3.48) 



has no solution. 
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^2,^1 



4 



Y 



Having this observation in mind our suggestion is to ignore this 
singularity. Our proposed strategy is the following: we ignore all such 
singularities that appear to be due to the possibility that two or more 
poles have the same position. Next we solve the problem and return to 
considering contributions to the scattering amplitude having the solu- 
tion both in the perturbative QCD region as well as in the saturation 
region. 

Finally, we consider only two contributions for the simplest 'fan' 



diagram: the first one, given by Eq. ( 3.35 ), is due to the exchange of 
two non-interacting Pomerons; and the second term of 3.35 which is the 
renormalisation of the amplitude of the Pomeron interaction with the 
target. Both these terms are the same as what we considered in the 
BFKL Pomeron calculus in zero transverse dimensions. 







Figure 7: The first enhanced diagram. 



3.4 The first enhanced diagram 

In this subsection we will calculate the simplest enhanced diagram shown 
in Fig. 1^, whose contribution has the form 



(3.49) 



du, du, d.3 du, .1 vl vl vl A(0.^i)A(0,^2)A(0,^3)A(0,.4) ^.(.o(y-r') ^ 

A(0, j/i)A(0, ;/4) 

g^(..)(y'-y") g^(.3)(y'-y") g^(.4)(y"-o) g*(fc^,g^,o,z/i) .g(fc;,g;,0,z.i) x 
+ 9'i/2 - 0, V2) g{-gi - qj2, - q^, 0, 1^2) X 

- 9;/2, - 0, z^a) g("fc4 + 9;/2 - ^3/2, 93, 0, 1^3) x 

dvi dv2 dv3 dvi vl vl vl vl A(0, V2) A(0, v^) e'^'y^^)^^ -^') x 
g.(..)(r'-y") ,^(.3)(y'-y") e^(..)(y"-o) <,*(fc^,^^, 0, z^O g(fc;, 9^,0,^.4) X 

+ 9i/2 - 9^/2, 0, V2) g{-kU - qj2 - ij2,qj2 + <[J2 - q^, 0, ^2) x 
3*(-fc'l - 9;/2, q^ - 0, v;) gi^k!^ + q^/2 - qj2 - ^J2, qj2 - x 

g*(fc4> 0,1^4) g(k^,q^-,^,Vi) 

where we used the fact that fc, — — ^4 + q'j2 — q 12, — —k'^ ^ ?t/2 and q„ = q^ — q'. ■ 

Considering, for simplicity, the case of the forward scattering (g^ — ' q^ = (^^ = q^ = Q and denoting ~ '^2 ^ 
as well as using the following notations: v\ = v^vt^ = v' and V2 = v-j, = va. u: representation this diagram 
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reduces to the following contribution 



A{Fig.§ - 



' diy Yl K^i) v''^ dvi dv' cPq' (fki (fk2 



X 5* (fc, 0,0,1.) 5(fci,0,0,j.) g*(-fc'i+g72,g',0,;/i).g*(-fc'i-g72,-g', 0,1/2) 
X gi~k^ + 972, g', 0, i^i) .g*(-fc2 - g72, 0, 1^2) 9* (k^, 0, 0, i^') g ik^,0, 0, 1^') 

1 r^'T 1 1 1 
X 7; — / awe -. r 7 r — - 



(3.51) 



The integral over u> can be rewritten as the sum over different contributions. We will show below that v — v' and 
the integral has the form 



1 

2TTi 



1 



1 



a+ioo 

duo e -. — r -. — r 



(3.52) 



Y 



1 



1 



Uj{v) — Uj{vi) — U0{V2) J ^{v) — ljj{vi) — Uj{v2) 



The first term leads to the renormalisation of the BFKL Pomeron intercept (the first term in the brackets) as 
well as of the vertex of the interaction of the BFKL Pomeron with the target (see Fig. ^) . The second term reduces 
to the exchange of two BFKL Pomcrons without any interaction between them. However, for 70 = ^ ~ ivo, given by 
Eq. ( |3.43 ), instead of Eq. (3.52), we have the pole of the third order which leads to the contribution 



2-Ki 



a + 200 



duj 



(3.53) 



We will not consider this contribution by the same reason that we used in calculating the 'fan' diagram since in this 
diagram two Ponierons with v and v' turn out to be outside of the saturation region while two Pomcrons with vi and 
V2 are located inside the saturation region where we cannot use the BFKL kernel to determine the values of their 
intercepts. Therefore, the first 'fan' and enhanced diagrams are very similar with respect to the integration over 
rapidity. However, in the case of the enhanced diagram we have an additional problem: the q' integration. Indeed, 
even if we assume that k > ko and g = 0, the integration over q' is restricted by the smallest of the two momenta k'l 
and /c2(see Appendix B for more details on q' integration). In other words we have the following region of integration: 



1. k > ki > k2 > kg. 

In this region the first term in Eq. ( [A-21 ) contributes for all four Pomerons and we have the product of two 5 
functions in the vertices, namely, S{1 + 7 — 71 — 72) <5(1 + 7' — 71 — 72). These leads to 7 = j'{h' = v') and the 
entire contribution looks the same as in the case of the first 'fan' diagram. It is worthwhile mentioning that 
for dipoles inside one BFKL Pomeron we have the same ordering in momenta. 

2. k > ki > kg but ki < /c2 while k2 > kg. In this kinematic region the first Pomeron with v and the fourth 



one with i/' have contributions which stem from the first term in Eq. (A-21). However, for two Pomerons with 



vi and 1^2 the second terms in Eq. (A-21) play the most important roles. Collecting all factors of fci one can 
see that the integration over ki has the form 



(3.54) 



dl exp ( (1 — 2 ^ ^ ^ ^^2) I 
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where I = ln{kf/kQ). The first term (1) in the bracket of Eq. (3.54) stemming from the integration over qi < kl, 
the second term (— ^ ~ ii^) originates from the contribution of the Pomeron with i^, the third as well as the 
fourth terms reflects the product of the factor (fc^)~2+*'^i from the upper vertex and factor (/c^)"^*"^^ from the 
lower vertex since qf — kf. The integration generates S{2 — j — ji — 72)- The integration over k2 has the same 
character as for the case of the first term and leads to (^(7' — 71 — 72). These 6 functions give 1 — 7 = 7'. 
The decomposition of Eq. (3.52) gives once more the terms shown in Fig. ^ and the pole of the third order 
of Eq. (3.53) since uj{l — 7) = ^(7) (see Eq. (3.7)). There is only one difference: in this pole of the third 
order, three Pomerons (with i/, i>i and V2) are inside of the saturation region. This fact does not influence our 
reasoning that we need not consider this kind of contribution. 



Finally, we can conclude that the first enhanced diagram has the same three contributions as the BFKL Pomeron 
calculus in zero transverse dimensions. 



3.5 Reduction of the emhanced diagrams to the system of non-interacting Pomerons 



In this section we will prove that an arbitrary enhanced diagram (see the diagram in Fig. can be reduced to the 
system of non- interacting Pomerons in the kinematic region: asY <C l/as. As far as uj integration is concerned 
this diagram can be written as 



a+ioo 



div 
2Td 



ujY 2 

e 7b 



^(7) 



m+l 



{IPI\lu) 



(3.55) 



where 75 is the vertex of Pomeron interaction with the dipole ^ S {1PI\(jj) is the one Pomeron irreducable diagram 
and can be written in the following form 



I](lP/|w) 



LU - cj(7i) - ^1^(72) 



1 



i>2 



Eu>2 ^ill) 



(3.56) 



to ■ 



In Eq. (3.55) and Eq. (3.56) 7 denotes the Pomeron - dipole vertex and F is the triple Pomeron vertex (see Eq. ( 3.25D ). 
The integral over uj in Eq. ( 3.55| ) can be taken closing contour over one Pomeron poles lo = 1-^(7) ( see the diagram 
in Fig. |l|-2) and closing contour over singularities of E (IP/jw) - the diagram in Fig. |^3. The first contribution has 
the form 

A{Fig. 0- 2) oc i- (y + O (1/as))" exp {iu{^)Y) E™ {\PI\u: = (3.57) 
'-' m! 

The contributions of the order of \/as stem from the differentiation of the factor E™ {1PI\ll}) with respect to w in 
Ec^. ( 3.55 ) in taking the pole of to + 1-order. All these contributions arc so small that they can be neglected, since 
Y > l/as (as^> 1). 

From Eq. ( |3.56| ) one can see that the first term in Ei {1PI\lS) is of the order of r^/w(7) (x Og while the others 
lead to smaller contribution. Therefore, the only term with to = remains that gives the exchange of one Pomeron 
(see Fig. |ll-4). 

The first singularity that Ei {1PI\lo) has is the pole uj — a;(7i) + 0^(72) . These poles are shown explicitly in 
Fig. ^3 and closing contour on these poles we obtain the contribution which is similar to Eq. ( 3.57 ), namely, 

A{Fig.%-^) (X —M{u;^Lo[^^)+oj{^2)/{Y + O (l/as))"' Ei (2P/|c^ = c^(7i + 1^(72)) (3.58) 
TO2! 



^In this section we use for this vertex the notation 73 instead 7 since we reserve 7 for the anomalous dimension. We hope that in 
other part of the paper we use 7 for the vertex will not lead to any misunderstanding. 
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1 

2) 



3) 



4) 



5) 



6) 



7) 



8) 



9) 



Figure 8: Reduction of tlie enlianced diagrams to a system of non-interecting Pomerons 



tj = tj(7i) +tj(72) 






C) 






u) = a;(7i) + Lo{^2) + ^^(73) 



Figure 9: Tlie vertices of two and tliree Pomerons interaction witli dipoles and E {2PI\uj = + cj(72)) 



where 777.2 is the number of the two Pomeron poles^; A'2 is the vertex of dipole-two Pomeron interaction (see Fig. ^ 
A) and Si (2P/|w) is the sum of two Pomeron irreducable diagrams. Since the first term in Si (2P/|w) gives the 
contribution which is proportional to (see Fig. ^B) we can neglect all contribution except with 7772 — reducing 
the diagram of Fig. ^5 to two Pomerons exchange of Fig. ||-8. To take the contribution from three Pomerons exchange 
in the diagram of Fig. ^3 we need to consider the diagram of Fig. ||-6 and close the contour in oj over the poles 
Lo = w(7i) +tj(72) +^(73). Repeating the same procedure and taking into account that for the vertex of dipole-three 
Pomerons we have the relation shown in Fig. we obtain the diagram of Fig. |^9. Continuing this procedure we 
see that we reduce the general enhanced diagram to the system of non-interacting Pomerons. 



^We assume that 71 and 72 are the same for all two Pomerons poles. Actually, only sum of 71 + 72 preserves. Strictly speaking we 
need to consider all t wo Po merons poles a bit different and take separately each of them. It is easy to see that in doing so we get the 
same result as in Eq. (3.58) after integration over relevant 7's 
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4. General solution for the simplified BFKL kernel 



4.1 The main idea 

The analysis of the first diagrams for the BFKL Pomeron calcuhis shows that they have the same kind of contributions 
as for the BFKL calculus in zero transverse dimensions and the new one that stems from the possibility that for 
some specific values of the anomalous dimension (70. fc) the intercept of n BFKL Pomerons in exchange is equal to 
the intercept of k BFKL Pomerons with k < n. However, for this specific contribution two or more Pomerons are 



located inside the saturation region where we cannot use Eq. (3.7) to calculate the intercept of the Pomeron. Inside 



the saturation region we can use Eq. ( 3.47 ) which has a rather general proof (see Ref. |^). In doing so we see that 
this specific QCD contribution leads to an energy suppressed scattering amplitude and can be neglected in the first 
approximation to the problem. 

Having this result in mind we can restrict ourselves by calculating only terms stemming from the decomposition 



of the type of Eq. (3.52). This decomposition leads to the exchange of non-interacting BFKL Pomerons and to the 
renormalisation of the vertex for the Pomeron target (projectile) interaction as well as to the renormalisation of the 
Pomeron intercept. In the kinematic region given by Eq. ( |l.l| ) we can neglect the renormalisation of the Pomeron 
intercept and find the renormalised vertex for the Pomeron target (projectile) interaction from the initial condition 
to our problem since it is closely related to the dipole -dipole (or target) interaction at low energies. 

Since we are dealing with the system of non-interacting Pomerons we can apply the procedure of the MPSI 
approximation (see Ref. [^) to take into account all the Pomeron loops. The strategy of the improved MPSI 
approach consists of three major steps 

1. To find the solution in the mean field approximation for arbitrary initial conditions. According to our analysis 
in this approximation the scattering amplitude has the following form 



N 



00 p 71 

{Y-Y',[jR{h,b,)]) = lld^k, i-^r Cn{k) P{k,h;b,\Y -Y') jR{h,h) (4.1) 

— 1 — 1 



where is an arbitrary function of k (recall that k is the conjugated variable to the size of a dipole). 

Comparing the exact solution with the initial condition to find the solution in terms of the renormalized 
vertices. It means that we need to solve the functional equation (in the general case) 

N{Y-Y' = 0,[jn{k,b)]-b) = N{[j{k,b)]) (4.2) 

where the function N should be given. In other words the scattering amplitude for low energies should be 
known either from calculations in QCD or from some phenomenology. This step is in full agreement with the 
parton model where the high energy interaction can be calculated through the low energy amplitude for the 
interaction of 'wee' partons that have to be given. 

To use the MPSI formula to take into account the Pomeron loops. This formula reads as follows 

A^"^^^^(fc,fco;fo,i^-0) = V TT d'k,d'k,d^,d%,^-^ 

^1 ^=\U '5^^^\hM)5^f{k„b,) 

N''^^ {y Y', N^"^^ (y', WR\k,,b,)]) l,a,^,<=)=„ 7^^ {h, k„b- b, - bj,) (4.3) 
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where 7^"* is the amphtude in the Born approximation. It should be stressed that Eq. ( }4.1| ) for 7V*^^^ 
guarantees that N^'^^psi (j^^ 5^ y _ q) does not depend on the vahie of Y' . Indeed, using the fact that 

J d^hdHjd\dHjPik,kf,h\Y-Y') ^^'^{h,k,,b^b,-h) P{ko,ky,b,\Y') = al P {k, ko;b\Y) (4.4) 



one can rewrite Eq. (4.3) in the form 



iV^^^^^(fc,fco;6,F-0) = V n d'k,d'k,d%d%^^ ^ 

N^^"" (Y-Y',[^^^\hM)]) N^'^^ {Y',[if{k„b,)]) l,a,^,(=^,^o {o^l P {k M-b\y)T 

= 1 - -P { - 4 P (fc, ^0; ^1^) ^ ^ } (^^'0 (^^'0 W-^.-^^. (4.5) 

where 7^,'^ = P (fc, fc,; - F') 7fl(fc., fe^) and 7^?^ = P {ko,h;h\Y') -fE{hM) 

At the moment we can carry out this program only for a simplified BFKL kernel since only for this kernel do we 
have the exact analytical solution to work with (see Q). The substantial number of numerical solutions Q does not 
help us since we know how to perform the renormalisation of the vertices only for the analytical solution. 

4.2 The simplified BFKL kernel and MFA solution 

The kernel for which we find the solution and implement our program has the following form |^ 

'"(^^ 2+^") = I ^ forr^Q^^i. (4-6) 

This kernel sums the leading log contributions in the region of low x. The first one is the usual LLA approximation 
of perturbative QCD in which the terms of the order of (as ln(r'^ A^))" are taken into account for ^ 1 

where as ln(r^ A^) 3> 1; and the second log approximation leads to the summation of the terms of the order of 
(as ln(r^ Qs))" in the kinematic region where » 1 and as ln(r^ Q^) > 1 ^. 

The LLA leads to the ordering of the dipole momenta. 

In the perturbative QCD region (r^ <C 1) we have, using the diagrams of Fig. ^ as an example, 

/c > fci > > fco (4.7) 

In the saturation region the ordering is the opposite, namely, 

/c < fci < < max( fco or Qs) (4.8) 

where Qs is the saturation scale. 



Eq. (4^) is not only much simpler than the exact kernel of Eq. (3/7) but calculating the Pomeron diagrams in 
LLA we can neglect the contribution due to the overlapping of the intercept for a different number of Pomerons in 
exchange. Indeed, the fact that a number of Pomerons are in the saturation region while other Pomerons are outside 



this region, means that we cannot keep the ordering of Eq. (4.7) or Eq. (4.8). Therefore, dealing with this kernel we 



enhance the arguments for neglecting the overlapping singularities. 
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Fortunately, the solution for the MFA has been found in Ref . in the entire kinematic region for arbitrary initial 
conditions. For completeness of our presentation, we will describe the main features of this solution inside of the 
saturation region in this section. The solution is written for the amplitude in the coordinate representation, namely 

'^^'^ .e't2Li2 N{k,b;Y) 



(2n) 



The solution in the MFA ( to the Balitsky-Kovchegov equation) can be written in the form 

N{z) = I - e'^(^) 

where the variable 

z = In {xl,QUY,b)) 
and function —({z) is determined by the following equation 



z 



= V2 f 
Jc, 



Co(b) + (exp(-C') - 1) 



(4.9) 

(4.10) 
(4.11) 

(4.12) 



The boundary conditions for the solution of Eq. ( l.lOp we should find from the solution at z < 0. We assume 
that to the right of the critical line Xi2Q'iiY,b) « 1 (2 = 0) the exchange of one BFKL Pomeron gives the main 
contribution. It has the form for our simplified BFKL kernel 



P{z,j) = 7 6 



7(6) 62 



(4.13) 



where ( in the spirit of the LLA approximation ) depends only on initial conditions at low energy and has a clear 
non-perturbative origin. Therefore the boundary conditions have the form 



N {z ^ 0+; Eq. (i.lO)) = iV(z = 0_) = 7 



lib); 



d In TV , 



dlnN. 



dz 



dz 



For small ( the solution to Eq. (4.12) has the form 



In- 



c 



z or C = Co(&)e^ 



(4.14) 
(4.15) 

(4.16) 
(4.17) 



Co(b) 2 

This equation says that we can satisfy Eq. (4.14) and Eq. ( 4.15| ) if 

Co{b) = 7(6) and 7(6) « 1 

Since 7(6) « a| Eq. ( [1.17 ) shows that Eq. ( 4.12 ) gives the solution to our problem. 

From Eq. (4.12) one can see that at small values of Co the integral has a logarithmic divergence. Therefore, we 
can rewrite Eq. ( 4.12| ) as follows 



z + 



2 ln7(6) = 2 ln(^7(fe)e^^) = ^ J 



dC 



Co(b) + (exp(-C') - 1) 



lnCo(&) 



(4.18) 



The r.h.s. of this equation does not depend on Co(^) = lib) (see Eq. ( 4.17 )). In Fig. |10| we plot ((z) for two cases : 
the exact expression of Eq. ( 4.12 ) and 



^/2 



dC 



a + (exp (-0 - 1) 



In (Co (6) /a) 



(4.19) 
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One can see that both functions coincide. This is an illustration that the integral of Eq. (4.12) has no further 



dependence on 7(6). Having this in mind we can claim that Eq. (4.18) leads to a general function, namely 



C(z) = ^ <f(7p(^,&)) 



(4.20) 



The argument of the function $ is the one Pomeron exchange (see Eq. ( 4.13| )) and expanding $ [j{b) ea ^ ) in a series 



with respect to this argument we are able to find the coefficient C„ in Eq. (4.1) 



Eq. (4.12) has also two simple analytical solutions for the function in the region of small ( ( see Eq. (4.16) 
and for large The last solution is 



C(z) = iln2(7(6)e^ 



(4.21) 



The exact solution is shown in Fig. [Til together with these two analytical solutions. One can see that these two 
analytical solutions give a good approximation for small and large values of 7(6) e 2 ^. 




4.3 Improved MPSI solution for Pomeron loops 



Using Eq. (4.1C) and Eq. (4.20) we can find the sum of the Pomeron loops by means of Eq. (4.5). We demonstrate the 
answer by solving analytically two extreme cases: (i) the behaviour of the scattering amplitude near to the saturation 
scale; and (ii) the asymptotic behaviour deep inside of the saturation region. 



For the first case ( — 7(6) exp (iz) and Eq. ( 4.1C ) together with Eq. ( |4.5| ) lead to the following answer 



N 



IMPSI 



(a;^i^2;r) = l-exp(-7(&)P(2;2;i?2.y)) = 1 - cxp (-7(6) e^^) 



(4.22) 



where x and R are the sizes of the projectile and target dipoles and z is defined by Eq. (4.11) where cx 1/ R 



For finding the behaviour of the scattering amplitude deeply inside of the saturation region we will use Eq. (2.18). 
In new variables: — h — In (^7p') ^^'^ + h = In (^7p'') j the contour C looks as it is shown in Fig. Since 
does not have singularities in the variables I and h we can replace the contour C by the contour Cr (see 



N 



MFA 
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Fig. p^ ). Therefore, we can rewrite Eq. (2.18) in the foUowing form 

1 

(27rz)- jc„ JCn 



N, 



1 



X Z 



MFA 



MFA I p (gll-h 



r 0, 



Using Eq. ( ff.lSD and Eq. (|4.2lD we rewrite 



jMFA 



1 - AT 



NFA 



(4.23) 



(4.24) 



where z = zi + Z2- 

One can see that the integration over li reduces to a very simple integral, namely 



1 



Cb. 



1 



00 



1 



e" / dlisin{2 7Tli) e '1 = 

2^71" 



erfi{'K) 



where erfi{z) = erf{iz)/i and erf{z) = E2{z) ~ /J^,exp (— t^) is the error function (see Refs. 
The remaining integral over / has the form 



MPSI 



(z) = 1 



(4.25) 
10). 

(4.26) 



erfi{'!T) 



dl \ \ — exp 



7"^^ g 2 ^ ^ 



r 



X sin (27r/) e » ' 



It is clear from Eq. (4.26) that the dominant contribution stems from I of the order of unity due to the exponential 
decrease at large I. Therefore, the behaviour of the integrand at large exp(iz) will dictate the resulting approach 



to 1 for the amplitude. At z » 1 the integrand gives 7(6) exp (- 
which leads to the following answer: 



V) behaviour. We can take the integral over I 



erfi(n) 

er/i(7r) 
2^7? vr 



o 

U |7(6)e5"-'| X sin(7r//2) e"^'' = 
T f ^/4 + i7r\ , ^./4 + i7T 



(4.27) 



7(6) 62- 



The solution of Eq. ( 1.27 ) shows the geometrical scaling behaviour and at high energies N — > 1 — exp(— ^z). 
In such a behaviour two features look unexpected: the geometrical scaling behaviour since the statistical physics 
motivated approach leads to a violation of this kind of behaviour [ pT[ ; and the fact that the asymptotic behaviour 
shows a very slow fall down in comparison with the expected N — 1 — exp (— C z^) This slow decrease 

shows that the Pomcron loops essentially change the behaviour of the amplitude in the saturation region. 

Our main difference to the statistical physics motivated approach is the fact that we took into account the 
impact parameter behaviour of the scattering amplitude. We observed that the integration in the Pomeron loops is 
dominated by the impact parameters of the order of the largest dipole size involved in two vertices assigned to the 
loop. We can face two different possibilities. In the first case the largest dipole belongs to the vertex , for which we 
assign our rapidity and for which the generating functional or other equation for our statistical system is written. 
This is the lower vertex of the Pomeron loop at the lower value of rapidity. In this case everything is in agreement 
with the statistical approach. In the second case the largest dipole belongs to the vertex at the top of the Pomeron 
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loop at higher rapidity value which is our future in terms of the statistical approach. The future behavior of this 
situation cannot be described by a Markov chain without taking into account the long range rapidity correlations. 
It should be stressed that the appearance of the overlapping singularities cannot be included in the statistically 
motivated formalism until the problem is reformulated in terms of the interaction of new effective particles related 
to these singularities. 







4.4 Self-consistency check 



Figure 12: The contours of integration in Eq. (1.23). 



As has been discussed we need to come back and consider 
the overlapping singularities that have been neglected. The 
check of self-consistency of our approximation means that 
we need to calculate all diagrams but replace the exchange 
of the BFKL Pomerons by the sum of all enhanced dia- 
grams (see Fig. 13) in the kinematic region where the BFKL Pomeron enters inside of the saturation region. For 
example, the simplest diagram of Fig. [t] has to be replaced by the diagram of Fig. |lj in the kinematic region that is 
related to the overlapping singularities. 

For the simplified kernel the value of 70 (see Fig. ||) is equal to 0.6 and this value is close to 7cr = 1/2. However, 
since 70 — 7cr does not depend on the value of rapidity, at high energies the solution that we need to check will be 
at z 3> 1. 

Indeed, the solution of Eq. (3.28) has the form of Eq. ( 4.13D only at small values of v. In the region of small v 
the expression of Eq. ( 3.2§| ) has the form 



p{kM\y) 



IV < 1 



cx 



di> exp 



exp 



1 



32 asY 



; asY 



(4.28) 



In the last line of Eq. ( 4.28|) we used the steepest decent method with the saddle point value of v 



^SP 



IQasY 



(4.29) 



Therefore, Eq. (4.13) for the Pomeron holds only for z <C 4 Vas^ and at i^sp <C 1. In this kinematic region 



vsp ^ 1- In other words, the region 7 
amplitude behaves as 



A 



7o corresponds to the large values of z where our scattering 
1 - (4.30) 



For the Pomeron near and inside of the saturation region we can rewrite the Pomeron diagrams in terms of the 
variables z and k introducing the Mellin transform for the Pomeron in the form 



1 i-a+ico Xz 

P(z) ^ — dX — — (4.31) 

27ri Ja^ioo A - A 



where the intercept Aq = -1. 
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Sum of all enhanced diasiams 




Figure 13: The graphic form for notation of wave bold Figure 14: The diagram that we need to replace the one 
line : the sum of all enhanced diagrams. in Fig. ^ to study the overlapping singularities. By bold 

line we denote the sum of diagrams of Fig. 



Eq. (4.30) has the form in A-representation: 



A = i - ^ 

A A + Ao 



(4.32) 



the coefficient C in Eq. ( 4.32| ) reflects the fact that we know the asymptotic behaviour inside the saturation region 
within the pre-exponential accuracy. 

In this representation the diagram of Fig. |lj will be proportional to 



-1 /•a+ioo 

A{F^g.^) (x — / dX 

J a~ioo 



1 



A- Ao 



C 



1 



A + Ao J A - Ao 



(4.33) 



instead of Eq. (||5|). 



One can see that we do not have the overlapping singularities. The main contribution stems from A = Aq and, 
since we have a double pole for this value of A the diagram describes the renormalisation of the Pomeron intercept 



which is proportional to q| and can be neglected in the kinematic region of Eq. (1.1) 



5. Conclusions 

The main results of this paper are: 



In the kinematic region of Eq. ( IT ) the BFKL Pomeron calculus after renormalisation of the Pomeron - target 
vertex, can be reduced to summing over the non-enhanced diagrams that describe the system of non-interacting 
Pomerons; 

The sum of enhanced diagrams for the system of non-interacting Pomerons can be calculated by means of the 
improved Mueller-Patel-Salam-Iancu approach, using an additional piece of information on dipole scattering 
amplitude at low energy (for example, using the Born Approximation of perturbative QCD for its determina- 
tion); 

For the simplified BFKL kernel given by Eq. ( |4.6D the analytical solution was found to have two unexpected 
features: the geometrical scaling behaviour is deeply in the saturation region and is slow to approach {N — > 
1 — C exp(— iz) instead of iV ^ 1 — Cexp(— z^/8) as it happens in the MFA formalism. 
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In this paper we confirm the result of Hatta and Mueller |Q that a new type of singularities appear in the 
BFKL Pomeron calculus: the overlapping singularities. They lead to a steeper asymptotic behaviour than the multi 
BFKL Pomeron exchanges give. It means that even the MFA formalism is doubtful without talking about summing 
over Pomeron loops. We argue that we can neglect these singularities since they stem from the kinematic region 
where two or more BFKL Pomerons are in the saturation region. We did the self-consistency check that such kind 
of singularities do not appear in the framework of the solution that we found. 

As has been mentioned several times in this paper, the geometrical scaling behaviour is in contradiction with 
the statistical physics analogy description and the Langevin equation based on this analogy[|l^. However, we need 
to recall that the statistical-like approach at the moment exists only as a QCD motivated model since the form of 
the noise term that has been used, is oversimplified. Indeed, we can prove that the BFKL Pomeron calculus can 
be described by the Langevin equation with some noise term but this noise is so complicated ( see |5^) that we do 
not know how to treat it. Instead of this the simplified model for the noise term has been used. This model, as 
far as we understand it, has its own difficulties: (i) it cannot describe the overlapping singularities that have been 
discussed here; (ii) it did not include the impact parameter dependence which is very important for any treatment 
of the Pomeron loops as it was shown here; and (iii) even for the toy model in which we neglect the dependence of 
the vertices on the dipole sizes, the statistical model predicts the constant cross section at high energies ^ while the 
exact solution (see Ref. ) leads to a decrease of the total cross section. 

We believe that the most vital problem, among these, is the correct impact parametre dependence. The calcu- 
lation of the enhanced diagrams show that the typical impact parameter in the Pomeron loops is of the order of the 
smallest size of the dipoles involved in the two vertices in the loops. It might or might not be the dipole from the 
vertex at rapidity Y for which we write the equation. If this dipole belongs to the second vertex, it is related to how 
we proceed from here and we do not think that it is possible to describe such a system in a statistical way, assuming 
only short range correlations in rapidity. The failure of the toy model description by the Langevin equation is , in our 
opinion, due to the fact that a statistical model cannot describe the renormalization of the intercept of the Pomeron 
and their applicability is restricted to the kinematic region of Eq. (^]^) . 

Since we consider the case of fixed coupling , even in deep inelastic scttering we have to take into acount all 
Pomeron loops. Indeed, as it was shown in Ref. 0| only running QCD coupling selects the 'fan' diagrams and leads 
to the mean field approximation for these processes. One cn see that the sum of Pomeron loopsat near to the 



saturation scale has a typical eikonal form (see Eq. (4.22)). This result supports the numerous saturation models 
that describes well the experimental data js^ . 

However, it seems more realistically, for a description of the deep inelastic processes first to solve the problem of 
summation of the Pomeron loops in the case of running QCD coupling. In this case we expect § the the mean field 
approximation will be able to describe the sunstantial part of kinematic region. More that ten years ago Braun and 
Levin suggested the procedure how to include the running QCD coupling in to th BKL equayion that describes 
the linear evolution in the region of high energy. This suggestion was based on the strong assumption (which by the 
way, were proved ) that the gluon reggeization being the essential part of the BFKL approach in the leading order 
will be preserved in the next-to-leading order as well. The formula that they obtained has a so called triumvirate 
structure, namely, for the basic process G(g) -I- G((f ) — > G{q — <f) the fixed coupling constand should be replaced by 



as 



Recently this conjucture has been proven (see Ref. [Q. This breakthrough allows us to include the running QCD 
coupling in the equation of the high density QCD. However, an additional problem arise that should be studied. 



^We got this information from S. Munier (lecture at GGI workshop on high density QCD) and from E. Naftali (private communication). 
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namely, such way of taking into account the running coupHng effect induces the contriutions of infrared renormalons 
(see paper of Levin in Rcf. | |5^ ) . Such a contribution cannot be treated perturbatively and can restrict the theoretical 
accuracy of our approach. Therefore, summing Pomeron loops with running QCD coupling is our next challenging 
problem. 
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Appendix A: Calculation of g (fc, q, n, i/) . 

We want to calculate the Fourier transform of the vertex function E defined in Eq. ( 3.13| ) as 



gih n, u) 



(A-1) 



±12 



with z — [x^ + 212)/-^ ~ 2o- "^^^ similar Fourier transform was already calculated in but the authors used a 
different normalization and we found it useful to calculate it using our notation. The integration over z was done 
by Navelet and Peschanski in by considering the solutions of the differential equations obeyed by the vertex 
functions E and consequent matching of the result to an approximate one obtained by Lipatov in . We adopt the 
notation introduced by Lipatov for the mixed representation of the vertex functions 



where 



On,u \p\ J 

2'^"'tt^ r(|n|/2 -iv+ l/2)r(|n|/2 + iv) 
\n\/2 - iv r(|n|/2 + w + l/2)r(|n|/2 - iv) 



and p = X12 in our notation introduced in Eq. (3.6). The expression for E''^'^{p) is given in | |43| and reads 
E'^'^i^p) ^ q ^'^-'^/^q ^'^+n/22-6^,.Y{l -iv+ |n|/2)r(l -iu- \n\/2) x 



(A-2) 



(A-3) 



(A-4) 



Thus in terms of E^''^(p) the Fourier transform in Eq. (A-1) is given by 
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g{k,q,n, u) 



(A-5) 



From Eq. ( A-4) and Eq. (A-5) one can see that the calculation of g{k, q, n, v) is reduced to evaluation of the following 
integral 



We introduce polar coordinates for complex variables 



k = \k\e 



iB 



and rewrite Eq. ( A-6 ) as 

1ti,ti'{q,k) 



pdp 



P Jo 



(A-6) 



(A-7) 



(A-8) 



For the sake of simplicity in Eq. (A-8) and further below we denote by p, q and k their respective absolute values. 
We use the series representation for the Bessel functions 



00 

•'-<=) = (I)" E 



(-!)• 



-^^ m\ T(hi + m + 1) 

7n—0 ^ ' 



(I) 



2m 



(A-9) 



to recast Eq. ( A-8 ) into 



00 f2-K 

dp 







^ipk cos((/>— ^) 



^iMC^-V-) ( 'IP 



(-1)^ 



m! r(/x + m+l) V 8 
(-1) 



(A-10) 



(qpy (--^j (qp\ 

V8/ ^^m'\r{n' + m' + l)\8 ) 



2m' 



-z2m' (0— ^) 



Changing variables (p — 9 

Jn{z) = 



1 

2^ 



and using the integral representation of the Bessel functions 

^ pTT-TT/2 



-in9+iz sin 6 



de 



271 



g-m(9-7r/2)+i2sin(e-7r/2)^^^ _ ^^2) 



(A-11) 



-7r-Tr/2 



2tt 



tt/2 



37r/2 



-in9-\-iz cos 9 



de 



2tt 



-inO-\-iz cos 



dO 



we rewrite Eq. ( A-10| ) as 

It,,t^'{q,k) = X] X] / '^P (^'^^ J-M-2m+/i'+2m(pfc) 



m— m'— ' 



jA''g-v'(S-i/') 



(f)' 



1 



to! r(/i + m + 1) 
1 



'gi2m(e--0) ^ 



(A-12) 



to'! r(/^' -I- to' -M) V 8 



QPY'"'~^2m'ie-i,) 



-29- 



Now we perform the integration over p using the formula 



I) 



(A-13) 



then the expression in Eq. ( |A-12| ) reads 



r?i=0 
00 



i! r(^ + TO + 1) V4fc 



1 



T{l/2- fi-2m) 



(A-14) 



- F 



7^ to'! r(^' + m' + 1) V4fc 



m'— 



(^■^y e-^2m'(e-V)r(i/2 + ;/ + 2TO') 



One should note the full factorization of two series in Eq. ( A-14| ) due to the integration over p. This allows us to 
sum them separately using 

,2m 



E 

rn—O 

and 



1 



_ cos(7r^) r(l/2 + fi) 
ml r(/i + m + 1) r(l/2 - p~ 2m) ~ n r(l + /i) 



[1/4 + p/2, 3/4 + /x/2; 1 + A*; 4 



E 



to! r(/i' + TO + 1) 



r(l/2 + p' + 2m) = + 2^^! [1/4 + A^'/S, 3/4 + p' /2- l + p';4z'] 

r(i + M ) 



(A-15) 



(A-16) 



With help of Eq. ( |A-15| ) and Eq. ( |A-16|) we obtain the final expression for I^,^' (q, k) as 



g \ cos(7r^) r(l/2 + 



4fc 



TT r(i + ^) 



2F1 [1/4 + p/2, 3/4 + p/2: l + p;4 ( A)%»2(e-^)] 



(^) r(l+7) 2Fr[l/4 + M72,3/4 + ^72;l + A^-;4 (A) (A-17) 

At this point we return to complex vector notation given in Eq. (A-7) and absorb the angles into complex vectors q, 
k, q and k 



^ I •-"(4)' 



1 



Til+p)T{l/2-p) 



2^^i[l/4 + A^/2, 3/4 + ^/2; 1 + /i; (Xy 



(A-18) 



^' f qy' r(l/2 + //) 



AkJ T{l + p' 



^ 2F1 [1/4 + p'/2, 3/4 + p'/2; l + p';[^)] 



In the first fine in Eq. ( [A- IS] ) we used identity r(l/2 + z)r(l/2 - z) = cos(z7r)/7r. Plugging Eq. (|A-18|) into Eq. (|A-5|) 
wc find the Fourier transform of the vertex functions E defined in Eq. (A-1) 



g{k, q, n,iy) ^ i"" q "'-''/^q «>^+«/22-6»>^r(l - ii' + |?i|/2)r(l - iiy - \n\/2) X 

— 27r 

[ln/2-iiy,-n/2-iu{q,k) — (— l)"I-ri/2+w,n/2+w (^j fc)] 



(A-19) 



with bn,i, given in Eq. ( A-3 ). In the high energy limit the main contribution to Pomeron Green function comes from 
n — 0. We are interested in the case of small momentum transferred along the Pomeron which implies 

|fc| » |g| (A-20) 

In the limit of Eq. ( |A-2C| ) for n = Eq. ( |A-19| ) reads 



gik,q,0,iy) = ^ \q\^"'2''^"'T^{l-iiy)cos{iiyTT)[J-,,,^,,{q,k)-I+,,^+,,iq,k)] 



(A-21) 



_22..|fe|-i+2.. r^(i/2-..)rM _ ^2-«.|^|-,-2..| I 



T'^{l + w)T{-w) 
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Appendix B: Triple Pomeron vertex in q and u representation. 



In this appendix we calculate the triple Pomeron vertex in q and v representation, namely, 
V^p{q[,q'^-v^,V2,vz) = J d?k[ g(fc'i,9'i,0,z.i) 9* {^^i + ^^.Q,^^ 3* - ^'i + ^^z'^i - 9^' 0' ^a) (B-1) 

where g is given by Eq. (A-21). 

Introducing new holomorphic variables 



t = 



Q'2 



2k[ + q!, 



and t 



Q'2 



2k[ + q'^ 



(B-2) 



and using Eq. (A-18) , we can rewrite Eq. (B-1) in the form 



Tap {0,q2,'^, 1^2,^3) = 4 ^ 



1/2— ii'i+ii'2+ii^3 



dtdt {{l-t){l-F))^/^+"'' Ciii^i) 



(B-3) 



X {Cl{v2)2Fi (1/4 + ^1^2/2, 3/4 + i;/2/2,l + it.2,t') 2F1 + iiy2/2,3/4 + 11^2/2,! + i,y2,P) 

- C*^{v2) [Utf'"-' 2F1 (1/4 + 11^2/2, 3/4 + 11^2/2, 1 + iiy2, t^) 2F1 (1/4 + IV2/2, 3/4 + 11^2/2, 1 + 1U2,?) } 
X { C*(i.3)2^^i (1/4 + ^^^3/2, 3/4 + iv^l2, 1 + zi^g, t') 2^1 (l/4 + iv^l2, 3/4 + 11/3/2, 1 + iv^, P) 

- C2*(z^3) {^tif 2F1 {I I A + iv^/2, 3/4 + ^I/3/2, 1 + iv^, P) 2F1 (l/4 + 11^3/2, 3/4 + ^^/3/2, 1 + ^l/3, } 



with 



CiM = —2 



2,, r^(i/2-z//)rM 
ii/" r2 (1/2 + ?;//) r(-ii/) 



C2H = — 2 



r2(i + izy) r(-ij/) 



(B-4) 



For simplicity we consider Fsp {q'i,q'2; vi, 1^2,1^3) in Eq. ( |B-3| ) for = 0. 



The integral of Eq. ( B-4 ) has three region of potential divergency: t ^ 0, t —>■ 1 and t 00. First, let us study 
the behaviour of the integrant at t — > 00. For doing this we use ( see formula 9.132(2) of Ref. |Q) the following 
expression 



2F1 (1/4 + iv2l2, 3/4 ^ 1V2/2, 1 + iV2,P 



r(l + zz.2)r(l/2) (l_ 
F2(3/4 + ii/2) V P 
F(l + zz/2)r(-l/2) 



2Fi ( 1/4 -I- W2/2, 1/4 + 1^2/2, ^ )(B-5) 



r2(i/4 + 11/2/2) 



4^^' 2 



2-Fi ( 3/4 + ^t/2/2, 3/4 + 11/2/2,^ 



From Eq. ( p3-5| ) we conclude that the integrant of Eq. (B-3) decreases at large t as for 1/2 = 1/3 = 0. It means 
that integrals over t and t are convergent as far as large t and t are concerned. The next suspicious region is i — > 
{t 0). Factor f)~i/2-H»i/i-«i/2-«i/3 ^^iq integrant indicates a possible divergence in this region. Rewriting the 
integral in a general form 



/>oo 

Tap (0,92; '^iif^2, 1^3) = / nd{tt){tt)' 
Jo 



l/2+ii'i— 11^2-21/3 



$(i,t) 



(B-6) 



d{tt) (tt) 



- 1 /2-\'il'i — 2f2 — i!^3 



{<i>(0) + (ci>(t,i)-$(0))} 



d(tF) (tt) 



-1/2 + 11/1 — ii/2— i!/3 



$(t,t") 
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we see that if $(t,t) — $(0) cx at t ^ we have a pole 1/(1/2 + ivi — iv2 — iv^)- However, since our function 
3>(t,t) has factors ti/2+"^i we need to be more carefuL In calculating first enhanced diagram (see section 3.4) V2 and 
are small and <i>(t, t) has the following form 



1>( 



(t,i) = Ci(i.i)^^ (l - {AtT)^'"'^) X (l - [Atlf^'-A 



(B-7) 



Using Eq. (B-7) we obtain that 



V2^0, V3- 



2 \ l/2-ii/i+ii/2+ii'3 



function without singularities in 1/2 and 1^3 

(B-8) 



(1/2 + ii/i - iv2 - iv^Y 

At — ^ 1 we could have a divergence which leads to a pole 1/(3/2 + wi). However, since Ci{vi) has a zero 
of the second order at ivi = —3/2 T^p (0, q2] I'l, 1^2, 1^3) has no such singularity. Therefore, Eq. (B-7) gives the triple 
Pomeron at small 1^2 and t/3. One can see that the pole in Eq. (B-8) replaces the 5- function which we obtained 
integrating first of q'2 ■ 

For arbitrary 1^2 and 1/3 the triple Pomeron vertex has a more complicated form, namely. 



92_ 



1/2 — 21/1 +11/2 +'it'3 



C*2{V2)CI{V3) 



Cl{y2)Cl{v3) 



Cl{v2)C*2{v:i) 



1/2 + ivi — iv2 — 11^3 1/2 + ivi — iv2 + iv^ 
C;{v2)CX{v3) 

^ ' regular function of V2 and V3 



{B-9) 



1/2 + ivi + iv2 — 11^3 1/2 + -I- iv2 + iv^ 

Integration over ((^ leads 8[y — v') where v (y2 in Eq. ( 3.49| )) corresponds to upper BFKL Pomeron in Fig. ^ 
and v' (1/4 in Eq. (3.49)) corresponds to lower BFKL Pomeron in Fig. 
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